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Abstract. Abundances of C and N are derived from features due to the CH G-band and to the UV CN band 
measured on high resolution {R <^ 40, 000) UVES spectra of more than 40 dwarfs and subgiants in NGC 6397, 
NGC 6752 and 47 Tuc. Oxygen abundances (or upper limits) are available for all stars in the sample. Isotopic 
ratios ^'^C/^^C were derived from the CH molecular band. This is the first determination of this ratio in unevolved 
dwarf stars in globular clusters. By enlarging the sample of subgiants in NGC 6397 studied in Gratton et al. (2001), 
we uncovered, for the first time, large variations in both Na and O also in this cluster. The origin of the chemical 
inhomogeneities must be searched for outside the stars under scrutiny. Our data indicate that in unevolved or 
slightly evolved stars in these clusters C abundances are low but not zero, also in stars with large N-enhancements 
and 0-depletions, and that the isotopic ratios ^^C/^^C are low, but never reach the equilibrium value of the 
CN-cycle. When coupled with the run of O and Na abundances, these findings possibly require that, in addition 
to CNO burning and p— captures, some triple— q process is also involved: previously evolved intermediate-mass 
AGB stars are then the most likely polluters. 
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1. INTRODUCTION 

After H and He, C, N, O are, together with Ne, the most 
abundant elements in the Universe. As such, they are key 
ingredients in a large number of astrophysical issues. Their 
abundances in metal-poor stars are tracers of the nucle- 
osynthetic sites that contributed to the different phases of 
galactic evolution. Moreover, they are important contrib- 
utors to the opacity in stellar interiors and act as catalysts 
in the CNO-cycle of H-burning. 

Stars of globular clusters (GCs) offer an ideal diagnos- 
tic in order to understand stellar evolution for low and 
intermediate stellar masses. However, since the pioneer- 
ing study of Osborn (1971) it is known that a spread 
in the light elements (C, N, O, but also among heavier 
species such as Na, Al and Mg) is present among cluster 
stars of similar evolutionary phase, unlike their analogs in 
the galactic field (e.g. Gratton et al. 2000 and references 
therein). Among these last, C and N abundances follow 
well defined evolutionary paths, with two episodes of mix- 
ing, the first one related to the first dredge-up and the sec- 
ond one after the red giant branch (RGB) bump, when the 
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molecular weight barrier created by the maximum inward 
penetration of the outer convective envelope is canceled 
by the outward expansion of the H-burning shell. The first 
episode was theoretically predicted by Iben (1964); while 
the second one is not present in canonical non-rotating 
models, it may be easily accomodated in models where 
some type of circulation is activated e.g. by core rota- 
tion (Sweigart & Gross 1978; Charbonnel 1994). Gratton 
et al. (2000) showed that among field stars no variations 
corresponding to these mixing episodes are observed for 
the remaining elements (namely, O and Na): again, this 
agrees with models, that do not allow deep enough mixing 
along the RGB. Clearly this pointed toward a peculiarity 
of globular cluster stars. 

As shown by Denisenkov & Denisenkova (1989), and 
later in a more quantitative way by Langer et al. (1993), 
the observed star-to-star scatter in cluster stars may be ex- 
plained by the CNO-cycle and the accompanying proton- 
capture reactions at high temperature. More uncertain is 
where these reactions occurred, whether in the observed 
star themselves, prior to an internal (extra- or enhanced-) 
very deep mixing episode, or elsewhere, perhaps in some 
form of H-burning at high temperature taking place e.g. in 
now extinct intermediate-mass AGB stars (IM-AGB), fol- 
lowed by ejection of polluting matter (see Gratton, Sneden 
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and Carretta 2004 for an updated review and references 
on the huge literature on this subject). 

Evidences from red giants are ambiguous, since both 
mixing (causing a decrease of [C/Fe] as a function of lu- 
minosity, see Bellman et al. 2001 and references therein) 
and pollution/accretion of processed matter (e.g. Yong ct 
al. 2003, Sneden et al. 2004) might be invoked to explain 
observations. 

Cleaner conclusions can be drawn from uncvolvcd 
or slightly evolved stars, where no mixing is expected 
and inner temperatures are not high enough to permit 
the p— capture reactions in the NcNa and AlMg cycles. 
However, up to a short time ago, only low-dispersion ob- 
servations of molecular bands of hydrides such as CH and 
NH or bi-mctallic molecules such as CN were available to 
study the chemical composition of faint GC turn-off (TO) 
dwarfs or subgiants (SGB). More important, no O indi- 
cator was accessible, since the atmospheric cutoff and the 
low throughput of existing spectrographs severely ham- 
pered the use of OH bands in the UV regions, and the 
remaining O features are only observable on high disper- 
sion spectra. 

In spite of these limitations, a number of studies (see 
e.g. Briley et al. 2004b and references therein) uncovered 
that large spreads in the CH and CN band strengths, an- 
ticorrelated with each other, do exist in unevolved stars in 
several GCs. The only explanation must necessarily rest 
on an event that polluted the material forming these stars, 
likely early in the cluster lifetime. 

Previous results from the present program (Gratton 
et al. 2001, Carretta et al. 2004) provided further 
strong evidences favouring primordial abundance varia- 
tions. Deriving the first reliable O abundances in cluster 
dwarfs, we found a clear anticorrelation between Na and 
O among TO and SGB stars in NGC 6752 and 47 Tuc. 
In NGC 6397, early results were not conclusive, but they 
were hampered by small number statistics. Similar anti- 
correlations were found also for Mg and Al. These facts 
require some non internal mechanism. 

However, after initial successes (see e.g. Ventura et al. 
2001), more recent models of metal-poor IM-AGB stars 
met serious problems in reproducing the 0-Na anticorre- 
lation and related phenomenology (Denissenkov & Herwig 
2003; Fenner et al. 2004; Herwig 2004), and unveiled that 
not only Hot Bottom Burning occurs (HBB), but also vig- 
orous H-burning at somewhat cooler temperatures dur- 
ing the intcrpulse phases. To overcome these problems, 
Denissenkov & Weiss (2004) recently proposed that the 
site for the p— capture reactions is the interior of RGB 
stars slightly more massive than those currently observed 
in globular clusters, and that they exchanged mass with 
the currently unevolved stars where the anomalous abun- 
dances are observed. 

In the present paper we complete the analysis of the 
spectra presented in Gratton et al. (2001) and Carretta et 
al. (2004) by including detailed abundances of C, N and 
O. From these abundances and from isotopic ^^C/^^C ra- 
tios, measured for the first time in such unevolved stars, we 



suggest that both triple-a captures in He-burning, to form 
fresh ^^C, and typical H-burning processing at high tem- 
peratures arc required to reproduce the observed pattern 
of abundances in these stars. If confirmed, this would ex- 
clude the possibility that mass-exchange with RGB stars 
might be responsible for the observed abundances. In the 
discussion, we will also comment on other possible short- 
comings of this hypothesis. 

Furthermore, we also present rcsiilts for Na and O 
abundances in a more extended sample of subgiants in 
NGC 6397, showing that large variations, anticorrelated 
with each other, in these two elements do exist also in 
this metal-poor cluster. 

2. OBSERVATIONS 

Details of observations are given in Gratton et al. 
(2001; Paper I) and Carretta et al. (2004). Briefly, spec- 
tra were acquired using the Ultraviolet- Visual Echelle 
Spectrograph (UVES) mounted at the ESO VLT-UT2 
within several runs (June and September 2000; August 
and October 2001, July 2002) of the ESO Large Program 
165.L-0263 (P.I. R. Gratton). On the whole, we have ob- 
servational material for 6 dwarfs and 9 subgiant stars in 
NGC 6397, 9 dwarfs and 9 subgiants in NGC 6752 and 
3 dwarfs and 9 subgiants in 47 Tuc. Relevant data for 
the observed stars are given in Gratton et al. (2001) and 
Carretta et al. (2004). Those for the additional subgiants 
in NGC 6397 are the same as given for the other subgiants 
in Gratton et al. (2001). 

Data were acquired using the dichroic beamsplitter 
#2. In the blue arm we used the CD2, centered at 420 nm, 
to cover both the CH G-band at ~ 4300 A and the CN UV 
system at 3880 A. The spectral coverage is about AA 
356-484 nm. The CD4, centered at 750 nm (covering AA 
555-946 nm), was adopted in the red arm. Observations in 
the run of Jime 2000 (mostly for NGC 6397) were made 
with a slightly different setup, resulting in a spectral cov- 
erage AA 338-465 nm in the blue and AA 517-891 nm in the 
red. Slit length was always 8 arcsec, while the slit width 
was mostly set at 1 arcsec (corresponding to a resolution of 
43000). In a few cases, according to the seeing conditions, 
this value was slightly modified downward or upward. 

In NGC 6397 and NGC 6752, typical exposure times 
were ^ 1 hour for subgiants, while in 47 Tuc we doubled 
this time. Each turn-off star was observed for a total of 
about 4 hours, split into several exposures. At A ~ 4300 
we reached a typical value of the S/N per pixel of ~ 30, 
increasing to ~ 70 for stars in NGC 6397. 

3. ATMOSPHERIC PARAMETERS 

The adopted values of the atmospheric parameters are dis- 
cussed in detail in Gratton et al. (2001) and Carretta et 
al. (2004). Here we only recall the main features of the 
analysis, that was also applied to the 6 newly observed 
SGB stars in NGC 6397. 
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We compared effective teniperatures from observed 
colours (both Johnson B — V and Stromgren b — y were 
used) with spectroscopic temperatures derived from fitting 
Balmer lines (namely Ha) . This approach was devised to 
derive precise values of the reddenings on the same scale 
for both cluster stars and field stars. These were used in 
the estimate of accurate distances to these clusters (see 
Gratton et al. 2003). 

Average values of temperatures were finally used for 
stars of NGC 6752 and NGC 6397. However, in 47 Tuc 
we found that the adoption of individual Tiff's for the 
subgiant stars provided best agreement with the values 
given by line excitation. Adopted values for 47 Tuc are 
given in Table 2 of Carretta et al. (2004). 

Values of the surface gravity were derived from the 
location of stars in the colour magnitude diagram; an age 
of 14 Gyr and corresponding masses were assumed. 

Estimates of the microturbulent velocity Vt for each 
star were derived, as usual, by eliminating trends of abun- 
dances with expected line strengths. Again, average values 
for each group of stars in similar evolutionary phases were 
adopted in NGC 6397 and NGC 6752. 

Finally, the overall model metallicities [A/H] were cho- 
sen as equal to the Fe abundances that best reproduce the 
measured equivalent widths {EW) , using the model atmo- 
spheres from the Kurucz (1995) grid with the overshooting 
option switched off. 

4. ANALYSIS 

4.1. Carbon and isotopic ratios 

Carbon abundances for the program stars were obtained 
from a comparison of observed and synthetic spectra in the 
wavelength region from 4300 A to 4340 A. This spectral 
region includes the band head of the (0-0), (1-1) and (2-2) 
bands of the A^A-X^H transitions of CH. We used newly 
derived line lists from Lucatello et al. (2003). Briefly, the 
starting line list was extracted from Kurucz 's database 
(Kurucz CD-ROM 23, 1995), including atomic species and 
molecular lines of C2, CN, CH, NH and OH. A few lines, 
missing in the Kurucz database, were added from the solar 
tables (Moore, Minnaert & Houtgast 1966); when uniden- 
tified, we arbitrarily attributed these lines to Fe I, with an 
excitation potential of EP=3.5 eV. 

The dissociation potential of CH has been determined 
with high accuracy at Dg — 3.464 eV (Brzozowoski et al. 
1976), while band oscillator strenghts were modified in or- 
der to reproduce the observed solar spectrum (Kurucz et 
al. 1984), using the solar carbon abundance of Anders and 
Grevesse (1989). We found that, in order to have a good 
match, a corrective factor of —0.3 dex in the \og gf val- 
ues of the electronic transitions and a shift of —0.05 A in 
wavelength were required, with respect to the values given 
by Kurucz. As found by many authors (see e.g. Grevesse 
& Sauval 1998), high excitation CH lines listed by Kurucz 
are missing in the spectra of the Sun and other stars, due 
to pre-dissociation. We omitted from our line list those 
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Fig. 1. The observed solar spectrum (Kurucz et al. 1984), 
shown as a continuous line, with overimposed the syn- 
thetic spectrum obtained using the adopted line list and 
the solar model (dotted line), in the G-band region at ~ 
4309-18 A. Notice that the spectral region shown in this 
Figure includes less than 1/4 of that used in our compar- 
ison with synthetic spectra. 



lines rising from levels with excitation potential over 1.5 
eV. 

The excellent match of the synthetic spectrum with 
the observed solar spectrum in part of the G-band region 
is shown in Figure ^ 

Using the appropriate atmospheric parameters, syn- 
thetic spectra in the spectral region 4300-4340 A were 
computed varying [C/Fe] in steps of 0.2 dex, in the range 
from 0.3 to —0.7 dex. A constant oxygen abundance 
[0/Fe] = was adopted in all these computations. The ex- 
act values affect only negligibly the derived C abundances, 
since for stars warmer than 4500 K the coupling of C 
and O is not relevant. 

After the synthesis computations, the generated spec- 
tra were convolved with Gaussians of appropriate FWHM 
to match the broadening mechanisms (in particular that 
due to the instrumental response) of the observed spectra. 
Carbon abundances were then derived from a set of 15-17 
CH features within the region under scrutiny, inspecting 
by eye all features, and computing an average value for 
each star. 

Results are summarized in Table ^ while Figure |21 
shows two examples of the synthetic spectrum fits to the 
observed CH features for a subgiant (star 478) and a turn- 
off star (star 1012) in 47 Tuc. The average rms deviations 
of the abundances from the individual observed features 
in [C/Fe] are 0.10-0.12 dex. 



4 



E. Carretta et al.; Abundance of C,N,0 in slightly evolved globular cluster stars 




Fig. 2. Left panel: spectrum synthesis of some features of 
CH band in a subgiant star of 47 Tuc. The heavy solid line 
is the observed spectrum, while dashed, solid and dotted 
lines are the synthetic spectra computed for three val- 
ues of the C abundances (listed on top of figure). Right 
panel: the same, for a dwarf star of 47 Tuc. Note that the 
synthetic spectra are now computed with different C val- 
ues. All synthetic spectra were convolved with a Gaussian 
to take into account the instrumental profile of observed 
spectra. 



Isotopic ratios ^"^C/^^C were estimated from spectrum 
synthesis in the two regions 4228-4240 A and 4360-4372 A 
containing various clean features of ^'^CH (see e.g. Sneden, 
Pilachowski and Vandenberg 1986; Gratton et al. 2000). 
Synthetic spectra were computed using the C abundances 
for each star derived from the G-band synthesis and the 
appropriate atmospheric parameters, and a range for the 
12q^13q values. The adopted isotopic ratios were derived 
as the averages from several features in both regions. 

4.2. Nitrogen 

Since the violet CN band at 4200 A is vanishingly weak 
in warm metal-poor stars (see Cannon et al. 1998, but see 
below the case of 47 Tuc) , we used the UV CN band and 
derived N abundances from a number of CN features in 
the wavelength range 3876-3890 A, where the bandhead of 
the Aw — bands of the UV system lies, again using line 
lists optimized by Lucatello et al. (2003). These lists use a 
CN dissociation potential of 7.66 eV from Engleman and 
Rouse (1975); a corrective factor of -0.3 dcx in the log g/ 
of the electronic transitions was applied also in this case, 
with respect to the values listed by Kurucz. In Figure|21the 
comparison between the observed solar spectrum (Kurucz 



Fig. 3. The observed solar spectrum (Kurucz et al. 1984), 
shown as a continuous line, with overimposed the syn- 
thetic spectrum obtained using the adopted line list and 
the solar model (dotted line), in the region at ^ 3880 A 
including the Av — bandheads of UV CN transition. 

et al. 1984) and the synthetic spectrum computed with 
the optimized line list is shown. 

Carbon abundances derived above from synthesis of 
the CH bands were adopted in the computation of syn- 
thetic spectra, relevant for individual stars, together with 
the appropriate atmospheric parameters (from Gratton et 
al. 2001 and Carretta et al. 2004). 

In the case of 47 Tuc, which is about 1.3 dex and 0.7 
dex more metal- rich than NGC 6397 and NGC 6752, re- 
spectively, we were able to use also the violet CN band 
strengths at ~ 4215 A in order to estimate the N abun- 
dances, at least in the subgiant stars. A procedure similar 
to that described above was used to compute synthetic 
spectra in the region from 4202 A to 4226 A, with the 
proper C abundance for each star. The N abundances re- 
sulted to be in very good agreement with those derived 
from the synthesis of the 3880 A region, so for the sub- 
giants in 47 Tuc the [N/Fe] values are those obtained as 
the average of N abundances in the two regions. 

No observations for the NH band were available for our 
program stars, apart from the very first run (June 2000), 
where the setup covered the region from 3376 to 3560 A, 
missed in the following observing runs. In this run stars 
in both NGC 6397 and NGC 6752, but not in 47 Tuc 
were observed. This choice of the setup was driven by the 
consideration that the expected S/N of the spectra of the 
(fainter) stars in NGC 6752 and 47 Tuc was so low that 
likely no meaningful abundances could be obtained. 

For these spectra of NGC 6397 stars, we then prepared 
a line list in the spectral range 3400-3410 A, where some 
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NH lines lie, using again the solar spectrum as a starting 
point; however, in order to obtain a good match, the log gf 
values of NH lines in Kurucz's list had to be lowered by 
about 0.5 dex. Results of the NH synthesis in stars of NGC 
6397 are given in the next Sect. 

4.3. Oxygen and Sodium 

Oxygen abundances in these warm stars were derived al- 
most exclusively from the permitted near-IR triplet at 
7771-75 A, as discussed at length in Gratton et al. (2001) 
and Carretta et al. (2004). Only for one subgiant in 47 
Tuc could we measure the forbidden [O I] lines. For the 
other stars, the very weak [O I] line was masked by much 
stronger telluric features, so that no reliable abundance 
could be derived. Final abundances and upper limits are 
given in Table ^ corrected for non-LTE effects as de- 
scribed in Gratton et al. (1999), from statistical equi- 
librium calculations based on empirically calibrated coUi- 
sional H I cross sections. The appropriate corrections were 
also applied to the Na abundances, derived from the strong 
doublet at 8183-94 A. 



5. RESULTS 

Derived abundances for C, N, O and isotopic ratios 
12^/13(3 foj. g^j^j,g 6397, NGC 6752 and 47 Tuc 

are listed in Tabled Carbon isotopic ratios could not be 
reliably derived for stars in NGC 6397 and dwarfs in NGC 
6752; only upper limits for C and N abundances were ob- 
tained for dwarf stars in NGC 6397, due to the weakness 
of the features and the low S/N ratio in the blue region 
of the spectra. 

Table n also lists, for an easier comparison of the rele- 
vant elements involved in H-burning at high temperatures, 
the abundances of Na taken from the previous papers of 
this series (Gratton et al. 2001, Carretta et al. 2004). For 
the 6 subgiants in NGC 6397 observed in July 2002, newly 
derived Na and O abundances are also shown in this Table, 
where stars are ordered according to increasing Na abun- 
dances. 

In the following, some features of the analysis of indi- 
vidual clusters are discussed. 

47 Tuc Being much more metal-rich than the other two 
clusters, 47 Tuc is the only one for which we were able 
to obtain meaningful lower limits of the isotopic ratios 
12q^13q (^^Q^j-f stars. The values found are listed in 
last column of Table ^ 

For the three dwarfs, the rather low quality of the 
spectra in the blue hampered a precise determination of 
the isotopic ratio. Hence, we choose to smooth somewhat 
the spectra, degrading the resolution to enhance the S/N. 
However, even in this case, the best result wc could secure 
is that the ratio ^^C/^^C is > 10 in these turn-off stars. 



It should be noticed that these are, to our knowledge, 
the first determinations, of the ^^C/^^C isotopic ratios in 
stars less evolved than the RGB-bump in GCs. 

NGC 6752 For dwarfs in NGC 6752, the quality of spectra 
does not allow a clearcut determination of the C abun- 
dances in the 0-poor dwarfs, which are very rich in N 
but with low C abundances. Hence, in these warm stars, 
the features of CH are rather weak, and we adopted the 
following procedure, in order to obtain a more reliable es- 
timate. 

The spectra of individual dwarf stars with low or 
not detected oxygen were summed up and this coadded 
spectrum was then used to derive an abundance of C 
through comparison with synthetic spectra. Our best es- 
timate is [C/Fe]= —0.2 ± 0.1. Analogously, the N abun- 
dances were then derived from the region 3876-3890 A 
using [C/Fe]= —0.2 dex and synthetic spectra computed 
with different values of the [N/Fe] ratio. Notice that ap- 
parently there are no N-poor dwarfs, and only one N-poor 
subgiant, in our sample. 

NGC 6397 In some stars of NGC 6397, acquired with the 
bluest setup during our first run, we were able to inves- 
tigate the NH molecular bands. The great advantage of 
using the hydrides bands in metal-poor clusters like NGC 
6397 is that bands of bi-metallic molecules like CN become 
vanishingly weak at low metallicity, due to their quadratic 
dependence on the metal abundance. Figure ^ shows the 
observed spectrum of the subgiant star 206810 as com- 
pared to five synthetic spectra computed with different 
[N/Fe] ratios. The lines of NH are clearly observed and 
the best match is obtained with the synthetic spectrum 
computed with [N/Fe]~ 1.3 dex, in very good agreement 
with the value that was derived from the CN bands. This 
supports our derivation of the N abundance for the sub- 
giants in this cluster. 

In Figure[Slthe same comparison is made with the aver- 
age spectrum obtained from the three dwarfs in NGC 6397 
having the best spectra in the UV, namely stars 202765, 
201432 and 1543. The resulting average spectrum was de- 
contaminated for a relevant (about 20% of the total value) 
contribution of scattered light due to the sky, not properly 
taken into account by our spectrum extraction procedure 
(note that these lines lie at the extreme UV edge of the 
observed spectrum). 

In this average spectrum we cannot firmly conclude 
that NH lines are actually observed: only some lines, but 
not all, are detected, and even these are very close to the 
noise level. The comparison with synthetic spectra in this 
region shows that a reasonable fit might be achieved at 
[N/Fc]< 1.5 dex: a value greater than 2.0 dex is clearly 
excluded. This upper limit is however more stringent that 
that derived from the CN lines. 
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Table 1. Abundances of C, N, O, Na, and isotopic ratios i^c/iSC in stars of 47 Tuc, NGC 6752 and NGC 6397 



Star 


fC/Fel 


[N/Fe] 


[O/Fe] 


[Na/Fe] 


^^C/'^C 






NGC 6397 - dwarfs 






1905 


<+0.50 


<2.0 


+0.24 


+0.09 




202765 


<+0.50 


<1.5 


+0.33 


+0.13 




201432 


<+0.50 


<1.5 


+0.21 


+0.15 




1543 


<+0.50 


<1.5 


+0.28 


+0.28 




1622 


<+0.50 


<2.0 


+0.23 


+0.35 








NGC 6397 


- subgiants 






706 


+0.10 


-0.5 


+0.54 


-0.48 




760 


+0.10 


-0.5 


+0.39 


-0.43 




777 


+0.15 


+0.2 


+0.39 


-0.25 




737 


+0.15 


+1.4 


<0.06 


+0.19 




793 


-0.10 


+1.2 


<0.06 


+0.23 




703 


+0.00 


+1.5 


<0.06 


+0.30 




206810 


—0.07 


+1.3 


<0.31 


+0.32 




729 


-0.10 


+1.4 


<0.06 


+0.48 




669 


+0.01 


+1.3 


+0.31 


+0.53 








NGC 6752 - dwarfs 






4428 


+0.09 


+1.1 


+0.33 


-0.29 




4383 


+0.12 


+1.2 


+0.57 


-0.18 




202316 


+0.12 


+1.5 


+0.27 


-0.06 




4341 


+0.21 


+1.4 


+0.20 


+0.20 




4661 


-0.20 


+1.4 


-0.35 


+0.29 




4458 


-0.20 


+1.5 


+0.02 


+0.31 




5048 


-0.20 


+1.5 


-0.30 


+0.37 




4907 


-0.20 


+1.5 


-0.25 


+0.58 




200613 


-0.20 


+1.7 




+0.62 








NGC 6752 


- subgiants 






1406 


-0.13 


+0.0 


+0.44 


+0.02 


9 


1665 


—0.28 


+1.0 




+0.10 


11 


1445 


-0.42 


+1.2 




+0.14 


5 


1400 


-0.25 


+1.0 


+0.38 


+0.20 


9 


1563 


-0.33 


+1.3 


+0.42 


+0.26 


5 


1461 


-0.32 


+1.3 




+0.28 


5 


202063 


-0.37 


+1.2 


+0.53 


+0.31 


3 


1460 


-0.49 


+1.4 


<0.29 


+0.42 


5 


1481 


-0.51 


+1.3 




+0.51 


4 






NGC 104 - dwarfs 






1081 


-0.13 


-0.50 


+0.57 


-0.34 


>10 


1012 


—0.10 


-0.30 


+0.48 


-0.14 


>10 


975 


-0.13 


-0.30 


+0.40 


+0.22 


>10 






NGC 104 


- subgiants 






482 


-0.16 


+0.10 


+0.61 


+0.06 


12 


206415 
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Values of [N/Fe] in the subgiants of 47 Tuc are the average of the abundances estimated from the synthesis of the 3883 A 
and the 4215 A regions. 
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Fig. 4. The observed spectrum of the subgiant 206810 in 
NGC 6397 (heavy soUd line) in the spectral region 3400- 
3410 A. The thin sohd hnes are synthetic spectra com- 
puted by using values of [N/Fe]= 1.0, 1.25, 1.50, 1.75 and 
2.0 from top to bottom, respectively. The NH lines are 
clearly observed. 

6. DISCUSSION 

6.1. The light elements: C and N 

Variations of the abundances of C and N in the examined 
stars, as derived from the features of the G-band and UV 
CN band, respectively, are shown in Figure Typical er- 
ror bars are also shown; they are conservative estimates, 
including both the scatter from the observed individual 
lines of CH and CN and the effect (almost negligible) of er- 
rors in the adopted atmospheric parameters (see Gratton 
et al. 2001 and Carretta et al. 2004 for the estimates of 
these uncertainties). 

In 47 Tuc the dwarfs are clustered around [C/Fe]^ 
—0.15 dex, [N/Fe]~ —0.4 dex, while the subgiants seem to 
be divided into two distinct groups, one with low N-high C 
and the other with low C-high N. Abundances of N and C 
seem to be anticorrelated in the other two clusters tt, even 
if for NGC 6397 we derived only upper limits for dwarfs, 
and the evidence of anticorrelation is somewhat weaker in 
the subgiants, with respect to stars in NGC 6752 and 47 
Tuc^ 

Even taking into account the small number statistics, 
it does not seem premature to conclude that in all 3 clus- 
ters there are a few subgiants with very low N abundances, 
well separated from high-N/low-C subgiants. The aver- 

^ Note that at [N/Fe]= -0.5 and [C/Fe]= +0.10 there are 
two subgiant stars in NGC 6397, where in Figure |S| only one 
point is displayed. 



Fig. 5. The same as in previous Fig. but for the average 
spectrum of 3 dwarfs in NGC 6397, namely stars 202765, 
201432 and 1543. 

age C abundance of the 3 low-N subgiants in 47 Tuc is 
[C/Fe]= —0.13 dex (rms = 0.03 dex, 3 stars), whereas 
the 6 subgiants with high N abundances have an av- 
erage of [C/Fe]= —0.34 {rms— 0.08). This difference in 
[C /Fe] is very similar to the one between the N-poor sub- 
giant in NGC 6752 ([C/Fe]= -0.13) and the average ob- 
tained from the other (N-rich) subgiants: [C/Fe]= —0.37 
{rms — 0.09, 8 stars). The spread in C abundances is 
smaller (about 0.15 dex) for subgiants in NGC 6397: in 
this case we obtain [C/Fe]= -1-0.12 dex {rms =0.12 dex, 
3 stars) and [C/Fe]= —0.02 dex {rms =0.10 dex, 6 stars) 
respectively for N-poor and N-rich stars. 

From these numbers we note what is immediately ap- 
parent in FigureEl while the spread in [N/Fe] is well above 
1 dex, in each cluster there is a rather small variation in C 
abundances. Since the C/N ratio in C-rich, N-poor stars is 
roughly solar (^ 0.6 dex), N in N-rich stars cannot be pro- 
duced only by transformation of C into N. Furthermore, 
even if carbon is a minority species in these stars, the 
residual C observed in N-rich stars is much more than 
that expected for material processed by CN-cycle at high 
temperature ([C/Fe] ^ — 0.8; see Langer et al. 1993) 

6.2. Comparison with previous works 

The anti-correlation between C and N abundances, al- 
ready known from low resolution spectra, is confirmed 
by our high dispersion spectra. To show this, we repro- 
duced in Figure |7| a similar plot shown by Briley et al. 
(2004a) with abundances of C and N for unevolved or 
slightly evolved stars in a number of clusters. While zero- 
point offsets are likely present between our data set and 
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Fig. 6. Run of the [C/Fe] ratio as a function of [N/Fe] 
for stars in 47 Tuc (red triangles), NGC 6752 (green cir- 
cles) and NGC 6397 (blue squares). Open symbols rep- 
resent dwarfs and filled symbols are subgiant stars, for 
all three clusters. Arrows represent upper limits in N, C 
abundances. Typical error bars are also shown. 



the [C/Fe] and [N/Fe] values derived by Briley et al. (as 
shown by mean ridge lines for 47 Tuc), the behaviour is 
essentially the same. In all clusters examined so far, varia- 
tions in C and N are anti-correlated with each other, with 
N showing large spreads, with respect to the more modest 
scatter in C abundances. Only among the SOB stars in 
M5 studied by Cohen et al. (2002) does the spread in C 
seem to equal the spread in N, and the most C-poor stars 
have a C depletion close to that expected by complete 
CN-cycling. 

Is this C-N anticorrelation tied to evolutionary phe- 
nomena occurring within the stars themselves or are we 
seeing the outcome of an already established nucleosynthe- 
sis implanted early in the material? To answer this ques- 
tion, we have to look into the evolutionary status of our 
program stars and seek for relationships with luminosity. 

In Figure|H|we plotted the [C/Fe] values of all program 
stars in the three clusters as a function of the absolute 
magnitude of the stars. Distance moduli are those cor- 
rected for the effect of binarity in Gratton et al. (2003). As 
a reference, we also plotted halo field stars from Gratton 
et al. (2000), in the same luminosity range; we only con- 
sidered field stars in the range —2.0 < [Fe/H] < —1, that 
closely matches the metallicity range of our three globu- 
lar clusters. This figure shows that the C abundance drops 
moderately (less than a factor of 2) both in field and clus- 
ter stars at the expected luminosity for the first dredge-up, 
in agreement with the standard stellar evolution models. 



Fig. 7. Run of the [C/Fe] ratio as a function of [N/Fe]. 
For our stars in 47 Tuc, NGC 6397 and NGC 6752 sym- 
bols are as in Fig. 6, with typical error bars shown. For 
literature data (all smaller symbols), filled yellow circles 
are SGB stars in M 5 from Cohen, Briley and Stetson 

(2002) , black crosses are main sequence turn-off stars in 
M 13 from Briley et al. (2004a), magenta empty triangles 
are M 71 turn-off stars from Briley and Cohen (2001) and 
magenta empty exploded stars are 47 Tuc MS stars from 
Briley et al. (2004b). 

Admittedly, the magnitude range is rather limited, but 
we can extend it by using literature data available for 
bright giants in the studied clusters. In FigureElwe added 
to our data also C measurements for red giants in NGC 
6397 (Briley et al. 1990), NGC 6752 (Suntzeff and Smith 
1991) and 47 Tuc (Brown et al. 1990). In all cases, the 
absolute magnitude scale is that defined by Gratton et al. 

(2003) . 

Despite the heterogeneity of data sources and methods 
used to obtain the C abundances (low-dispersion spectro- 
scopic measurements of the G-band in stars of NGC 6397 
by Briley et al. 1990; infrared spectra of the first overtone 
CO bands for stars studied by Suntzeff and Smith I99I 
in NGC 6752; moderately high resolution spectra of the 
CN red system and of the G-band for stars in 47 Tuc by 
Brown et al. 1990; synthesis of high resolution spectra of 
the UV CN system and of the G-band in our program 
stars), and offsets between the bright giants and the un- 
evolved stars might be present, the overall pattern among 
cluster stars seems to follow rather well the one defined by 
"undisturbed" field stars. Namely, the second step-like de- 
crease in [C/Fe] ratios around the red giant branch bump 
(A/y ~ 0.5) seems to be present also for cluster stars. This 
is explained (e.g. Sweigart and Mengel 1979; Charbonnel 
1995; Gratton et al. 2000) as the onset of a second mix- 
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Fig. 8. Run of the [C/Fe] ratio as a function of absolute 
magnitude for program stars in NGC 6397, NGC 6752 and 
47 Tuc. Symbols are as in Fig. 6 . Small (magenta) crosses 
are the field stars of the Gratton et al. (2000) sample with 
metallicity in the range —2.0 < [Fe/H] < —1. 

ing episode during the red giant evolution of a population 
II star, once the molecular weight barrier established by 
the retreating convective envelope is wiped out by the ad- 
vancing shell of H-burning. From now on, CN-processed 
material is able to reach the surface layer, where a further 
decrease of C is visible, as shown by Figure]^ 

We conclude that in spite of a larger (intrinsic) scatter, 
the same mixing episodes observed among field stars can 
be traced also among cluster stars. 

This result is not new. Very recently, Smith and 
Martell (2003) used the same field sample by Gratton et 
al. (2000) and literature data for C abundances in red gi- 
ants in M 92, NGC 6397, M 3 and M 13 to show that 
the rate of decline of [C/Fe] on the RGB as a function 
of My is very similar between cluster and halo field gi- 
ants. Our study, however, has the advantage of sampling 
regions along the RGB well below the so called bump in 
the RGB luminosity function, where standard theories for 
extra-mixing (see e.g. Sweigart and Mengel 1979) fix the 
threshold in magnitude for the onset of additional mixing. 
In globular clusters the chemical anomalies can be traced 
down to very faint magnitudes and we clearly detect a 
steady increase in the average C abundance going from 
red giants to subgiants and to dwarf stars. 

FigureElalso shows another well known feature: cluster 
stars seems to reach more extreme C depletions than those 
experienced by field analogs, as clearly indicated by red 
giants in NGC 6752 and, partly, by giants in NGC 6397. 

The analogous plot for N abundances, in Figure ^1 
also reveals an extreme behaviour of cluster stars. Spreads 



Fig. 9. Run of the [C/Fe] ratio as a function of absolute 
magnitude for program stars in NGC 6397, NGC 6752 and 
47 Tuc. Symbols are as in Fig. 6. Small (magenta) crosses 
are the field stars of the Gratton et al. (2000) sample with 
metallicity in the range —2.0 < [Fe/H] < —1. Open (blue) 
diamonds are the bright giants in NGC 6397 from Briley et 
al. (1990); open (green) diamonds are red giants in NGC 
6752 from Suntzeff and Smith (1991); open (red) diamonds 
are red giants in 47 Tuc from Brown et al. (1990). 



in [N/Fe] are very large in cluster dwarfs and subgiants 
with respect to field stars of similar evolutionary status. 
Moreover, when coupled with literature data (from Brown 
et al. 1990), a hint for increasing [N/Fe] at increasing lu- 
minosity seems to appear for stars in 47 Tuc. On the other 
hand, no clear indication of such increase in NGC 6752 is 
present. 

Finally, Figure [TTl shows the isotopic ratios ^^C/^'^C 
measured in program stars as compared to the field 
database of Gratton et al. (2000), as well as the literature 
data available for these clusters. Also in this case, clus- 
ter stars show extremely low values of the isotopic ratios, 
when compared to the field stars of similar magnitude. 
However, the interesting feature here is that the ^^C/^'^C 
values are not at the CN-cycle equilibrium value, not even 
for stars that are extremely N-rich. For these stars one 
would expect a value of ~ 3, at odds with our findings. 
Even considering literature data for bright giants, their 
isotopic ratios seem to be somewhat lower than those in 
field stars but always at a level slightly higher than the 
equilibrium value. 
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Fig. 10. Run of the [N/Fe] ratio as a function of absolute 
magnitude for program stars in NGC 6397, NGC 6752 and 
47 Tuc, field stars (Gratton et al. 2000) and cluster red 
giants from literature (Brown et al. 1990). Symbols and 
color codes are as in the previous Figure. 

6.2.1. First conclusions from light elements 

In summary, the analyses of C and N and the relative 
abundances of the carbon isotopes in slightly evolved glob- 
ular cluster stars show that: 

— there are moderately small variations in C abundances, 
anticorrelated to (large) variations in N abundances. 
However, except for a few stars in M5, the C is gener- 
ally not as low as expected for CN-cycling material. In 
no dwarf or early subgiant do we find C abundances 
as low as those observed in (all) highly evolved RGB 
stars, i.e. stars brighter than the RGB-bump. 

— N shows a large spread and, with the cautionary warn- 
ing of a rather limited range in sampled magnitudes, 
does not seem to have a particular increase at the 
first dredge-up position. Note that a large fraction 
of dwarfs and early subgiants have N abundances as 
high as those observed in highly evolved RGB stars 
(Figure Cni). 

— the isotopic ratios ^^C/^'^C are low but do not reach 
the equilibrium value of 3. 

From these observations, it is already clear that no un- 
evolved star has the same surface composition as highly 
RGB stars, as required by some recent scenarios advo- 
cated to explain the 0-Na anticorrelation (Denissenkov 
and Weiss 2004). Rather, the observed pattern resembles 
much more that predicted by the evolutionary models of 
the most massive intermediate-mass, low metallicity AGB 
stars (e.g., Ventura et al. 2002, 2004). According to these 



Fig. 11. Run of the isotopic ratio ^^C/^^C as a function 
of Mv for program stars in NGC 6397, NGC 6752 and 
47 Tuc. Symbols are as in the previous Figures. Small 
(magenta) crosses and upper limits are the field stars of 
the Gratton et al. (2000) sample with metallicity in the 
range —2.0 < [Fe/H] < —1. Open (green) diamonds are 
red giants in NGC 6752 from Suntzeff and Smith (1991); 
open (red) diamonds are red giants in 47 Tuc from Brown 
et al. (1990) and Shetrone (2003). 



models, these stars have complete CNO cycling at the base 
of the convective envelope and the coupling of nucleosyn- 
thesis plus release of processed matter into the protoclus- 
ter clouds or onto second-generation stars is able to pro- 
duce a surface composition where C is depleted, but not 
to values as low as those expected from CN-cycling, be- 
cause some fresh ^^C produced by triple— a reactions is 
dredged up to the surface, N is enhanced and the carbon 
isotopic ratio approaches the equilibrium value, due to the 
large enhancement of ^^C. Apparently, this is almost ex- 
actly the chemical composition of the unevolved cluster 
stars under scrutiny. 

At this point of our discussion, however, this asser- 
tion is not yet proven, because dilution with material 
not contaminated by CNO burning could be invoked to 
explain the observed trends for C and N abundances. 
Observations of heavier elements involved in high tem- 
perature p— capture reactions may give a deeper insight. 



6.3. Oxygen and Sodium 

Looking now at heavier elements (O and Na), we proceed 
along a path of stronger Coulomb barriers. The temper- 
atures involved arc much higher and we are considering 
deeper regions in the H-burning shell. 



E. Carretta et al.: Abundance of C,N,0 in slightly evolved globular cluster stars 



11 



The well known Na-0 anticorrelation (see Gratton, 
Sneden and Carretta 2004 for a recent review) is sum- 
marized in Figure 1121 for our program clusters. In this 
Figure we also added the available literature data, which 
is mostly for bright red giant stars, even if, apart from a 
few cases, no systematic studies have been performed for 
these 3 clusters, often used as calibrators. Abundances 
of Na and O in NGC 6397 include 2 stars studied by 
Norris and Da Costa (1995) and 2 stars from Castilho et 
al. (2000). For 47 Tuc, additional data are from Norris and 
Da Costa (1995) and Carretta (1994). For NGC 6752, we 
added 6 stars from Norris and Da Costa (1995), 4 stars 
from Carretta (1994) and the bright red giants studied 
by Yong et al. (2003). Individual values of Na and O for 
the bump stars in NGC 6752 analyzed by Grundahl et al. 
(2002) have not been published anywhere, therefore they 
are not used. Whenever possible, as in the Yong et al. 
sample, we started from original logn(O) and logn(iVa) 
values, bringing them onto a homogeneous scale by using 
the average [Fe/H] values and the solar values adopted in 
the present study. 

Note that values for our dwarfs and subgiants do in- 
clude corrections for departures from LTE. For literature 
data, this was possible only for stars analyzed in Carretta 
(1994). However, since O abundances are usually derived 
in red giants from the forbidden [O I] doublet, these correc- 
tions are negligible. NLTE corrections for Na abundance 
might be of more concern in giants, depending on what 
lines were used in the various analyses, but the overall ap- 
pearence of Figure IT^ shows that if there are some offsets, 
they are rather small. 

The 0-Na anticorrelation is in fact very well defined for 
all clusters; there seem not to be large differences among 
the different clusters over the metallicity range sampled, 
nor among stars of different evolutionary status within a 
given cluster. For the first time, the existence of a Na-0 
anticorrelation also among stars in NGC 6397 is clearly 
shown. This Figure shows among slightly evolved cluster 
stars the same trends that have previously been observed 
among the red giant stars of several globular clusters. 

The difference with respect to field stars is striking. 
With their highly homogeneous sample, Gratton et al. 
(2000) convincingly showed that field stars have constant 
Na and O abundances, not anticorrelated with each other. 
The obvious implication is that in cluster stars there 
is something else able to alter simultaneously the abun- 
dances of these two elements. 

To get a deeper insight, we show in Figure 1131 and 
FigureElthe derived abundances of C and N, respectively, 
for program stars as a function of the [Na/Fe] ratio. 

In all three clusters. Figure El and Figure El clearly 
show a trend for C and N abundances to decrease 
and increase respectively with the increase of Na abun- 
dances. In particular, the C-Na anticorrelation closely 
mimics the well-known 0-Na anticorrelation, summarized 
in Figure El 

Furthermore, while turn-off stars show a large range 
in Na abundances (at almost constant C), carbon abun- 
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Fig. 12. Run of the [Na/Fe] ratio as a function of [0/Fe], 
for stars in 47 Tuc, NGC 6752 and NGC 6397. Symbols 
for our program stars are as in Fig. 6. Literature data are 
as follows: (green) diamonds with crosses inside are bright 
red giants from the extensive study by Yong et al. (2003), 
open (blue, green and red) diamonds are stars of NGC 
6397, NGC 6752 and 47 Tuc, respectively, from Norris 
and Da Costa (1995), Carretta (1994) and Castilho et al. 
(2000), as described in the text. 
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Fig. 13. Run of the [C/Fe] ratio as a function of [Na/Fe], 
for stars in 47 Tuc, NGC 6397 and NGC 6752. Symbols 
are as in Fig. 6 . 
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Fig. 14. Run of the [N/Fe] ratio as a function of [Na/Fe], 
for stars in 47 Tuc, NGC 6752 and NGC 6397. Symbols 
are as in Fig. 6 . 
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Fig. 15. Run of the [N/Fe] ratio as a function of [0/Fe], 
for stars in 47 Tuc, NGC 6752 and NGC 6397. Symbols 
are as in Fig. 6 . 



dances are anticorrelated with Na abundances for SGB 
stars. On the other hand, N abundance correlates well with 
sodium among subgiant stars, even if the anticorrelation is 
less evident among TO stars. Finally, O is anticorrelated 
with N, as shown in Figure [T5l 



The overall distribution of light elements seems to 
point out that processes of proton-capture are at work. 
In the atmospheres of the stars studied we are just seeing 
the products of these reactions. In this case, the line of 
thought is the same as in Gratton et al. (2001): turn-off 
stars do not reach the temperature regime where the ON 
and NeNa cycles required to produce the Na-0 anticorre- 
lation are active, and moreover these stars have convective 
envelopes that are too small to have efficiently mixed the 
ashes of these nuclear processes up to the surface. The 
same conclusion holds also for subgiants. 

The bottom line is that what we are seeing are the 
by-products of nuclear burning and dredge-up in other 
stars, that are now not observable, but that have returned 
their elements to the intracluster medium or directly to 
the surface of the presently observed stars. 

6.4. CNO arithmetic 

Having now also O abundances at hand for program stars, 
it is possible to test in another way if the observed pattern 
of C, N, O abundances can be explained by the CNO-cycle 
alone. In fact, in this hypothesis it is only the relative 
content of C, N and O that may change, as a consequence 
of different reaction rates; their sum has to be constant. 

In Figure ITHI the C abundances are plotted as a func- 
tion of the sum C+N for our program stars and, as a 
reference, for stars with abundances from low dispersion 
studies, as in Figure |7| 

From this Fig. one has the impression that the sum 
C-l-N increases with decreasing C abundance; this is con- 
firmed by computing Kendall's r, which implies that this 
anti-correlation is highly significative, at the 99.3% level. 
Hence the observed pattern cannot be due simply to a 
transformation of C into N by the incomplete CN-cycle of 
H-burning. 

The possible alternatives are then either that the ON- 
cycle is also involved, adding N freshly produced from O, 
or that some variable amount of already existing N is su- 
perimposed on the effects of C-^N reconversion. In this 
regard, clearcut evidence is provided by Figure [T7I where 
the run of [N/Fe] as a function of the total sum C-t-N-|-0 
is shown. Over a spread in N of almost 2 dex, the sum 
remains almost (but not exactly, see below) constant; this 
by itself implies that the complete CNO-cycle has been 
at work to produce the observed pattern. Once more, nei- 
ther subgiants nor, in particular, unevolved turn-off stars 
are able to forge elements (such as Na) that require high- 
temperature proton-capture reactions. Moreover, they are 
unable to mix in their convective envelopes such products 
to the atmospheric layers; hence, we are forced to conclude 
that the CNO cycle was at work in stars other than those 
presently observed. 

The best candidates at hand are the intermediate- mass 
AGB stars. 

Another class of possible polluters was recently sug- 
gested by Denissenkov and Weiss (2004). According to 
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Fig. 16. Run of the [C/Fe] ratio as a function of 
[(C+N)/Fe]. For our stars in 47 Tuc, NGC 6397 and NGC 
6752 symbols are as in Fig. 6. For literature data (all 
smaller symbols), filled yellow circles are SGB stars in M 5 
from Cohen et al. (2002), black crosses are main sequence 
turn-off stars in M f 3 from Briley et al. (2004a), magenta 
empty triangles are M 71 turn-off stars from Briley and 
Cohen (2001) and magenta empty exploded stars are 47 
Tuc MS stars from Briley et al. (2004b). 




[0/N] 

Fig. 18. Run of the [C/N] ratio as a function of [0/N] for 
our stars in 47 Tuc, NGC 6397 and NGC 6752; symbols are 
as in Fig. 6. Superimposed on the data are the three mod- 
els outlined in the text: a simple dilution with material 
processed through the complete CNO cycle (solid, black 
line), contamination from N-poor RGB stars (dashed blue 
line) with composition typical of field RGB stars, and 
contamination from N-rich upper-RGB stars experiencing 
very deep mixing (dotted magenta line). 




Fig. 17. Run of the [N/Fe] ratio as a function of 
[(C-hN-hO)/Fe] for our stars in 47 Tuc, NGC 6397 and 
NGC 6752; symbols are as in Fig. 6. 



their computations, as well as those previously reported 
by Denissenkov and Herwig (2003), nucleosynthesis in IM- 
AGB stars with strong 0-depletion is not accompanied 
by large Na production (hence, the matter is not Na- 
enhanced as required by the Na-0 anticorrelation); in- 
stead, strong Mg depletions are expected, and this has 
never been observed in globular cluster stars. Similar re- 
sults has been recently obtained by Herwig (2004) and 
Fenner et al. (2004). As a way out, Denissenkov and Weiss 
(2004) suggested that pecuhar CNO abundances, as ob- 
served in unevolved cluster stars, might be a result of 
the H-burning shell in upper RGB stars of mass slightly 
larger than those presently observed in GCs, provided that 
they have experienced some degree of extra-mixing (see 
Denissenkov and Herwig 2003), followed by mass transfer 
onto less evolved stars. 



6.4.1. A simple dilution model 

So far, it has been shown (i) that globular cluster stars 
exhibit the same mixing episodes observed for field stars; 
(ii) that no slightly evolved cluster star has an abundance 
pattern the same as that observed among stars close to 
tip of the red giant branch; and (iii) that there may be 
an excess in the sum of C-l-N-l-0 in N-rich stars, that can 
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possibly be attributed to some ^^C in excess with respect 
to the predictions of the complete CNO cycle. 

To test how various schemes work to explain these 
observational facts, we will consider here simple models 
for the dilution of the products of CNO burning at vari- 
ous temperatures, and will compare their predictions with 
the run of the [C/N] ratio against [0/N] ratio shown in 
Figure El as well as with the other pieces of information 
we have gathered. 

When constructing our dilution models, we started 
noticing that in the complete CNO-cycle at high temper- 
ature (40 X 10^ K), at equilibrium the O abundance is 
decreased much more than the C abundance (factors of 
about 50 and 6, respectively: Langer et al. 1993). This is 
quite different from the case of the incomplete CN-cycle 
at low temperature (~ 10^ K), where the C abundance is 
decreased by a factor of 6, as before, but the O abundance 
is not modified. 

Testing the scenario of pollution by the complete 
CNO-cycle 

We started by considering the standard scenario of pol- 
lution by the complete CNO-cycle. Let a be the fraction 
of gas forming uncvolved stars in globular clusters that 
has been burnt through the complete CNO-cycle and let 
us assume the remaining gas to be of primordial com- 
position. It follows that the abundance of C and O in 
the mixed gas with respect to the initial composition is: 
f = a ■ b{C, O) + {1 — a), where b{C, O) is the abundance 
of C and O at the equilibrium of the CNO-cycle (0.17 and 
0.02 for C and O, respectively). The N abundance can be 
derived by the constraint that the total CNO abundance 
is constant. 

If we assume that initially [C/N]=0.5 and [O/N]==1.0 
(which is approximatively the composition of N-poor stars 
in the three clusters: see Figure [TH|l . we may then predict 
the values of [C/N] and [0/N] for different fractions a of 
the gas consumed in the complete CNO-cycle, with a — 
being the original composition. This trend is shown as a 
solid line in Figure El This line reproduces fairly well 
the location of observational points in Figure El albeit 
it predicts somewhat too low C abundances for N-rich 
dwarfs. Dilution factors a adequate to reproduce N-rich 
stars are 0.5 < a < 0.8. Such values also allow reproducing 
the isotopic ratio ^^C/^'^C~ 7 observed in unevolved stars. 
Moreover, with these factors we might quite easily obtain a 
Li abundance roughly similar to the original one, provided 
that the diluting material be Li-rich (as expected from 
some of the IM-AGB stars: Ventura et al. 2001). 

On the other hand, with such large dilution factors, the 
models by Langer et al. (1993) would predict too much 
Na, by more than a factor of about 10; this is because 
in these models Na is produced also by ^°Ne. In order 
to reproduce the observations we would need that Na be 
forged only from ^^Ne. 

Testing the scenario of polluting RGB stars 



Let us now use a similar model to test the scenario 
envisioned by Denissenkov and Weiss (2004) of pollution 
by RGB stars. In this case, mixing occurs in a fraction 
of upper-RGB stars and afterward a transfer of material 
onto the dwarfs occurs. 

Let us assume that the upper-RGB stars might 
have either one of the following 2 compositions: (i) a 
chemical composition typical of field upper-RGB stars 
(N-poor, i.e. N only coming from incomplete CN- 
cycle) or (ii) a composition from very deep mixing, 
where complete-CNO cycle and Na-enrichment are in- 
volved. For these stars we will use the most extreme 
case observed (i.e. N-rich). Note that in both groups 
of stars all Li is destroyed. The starting composi- 
tions (logn(C)/(N)/(0)/(Na)/(C-f N-KO)) assumed are 
then: 8.60/8.00/8.90/6.30/9.11, 8.60/7.50/9.30/5.80/9.38, 
7.50/9.33/8.40/6.80/9.38 and 8.00/8.50/9.30/5.80/9.38 
for the solar, original, N-rich and N-poor cases, respec- 
tively. 

By varying the dilution factor a, we obtain the abun- 
dance pattern for the two cases original -|-a x (N-rich) and 
original-l-ax (N-poor). Results are overplotted as a dotted 
(magenta) line and a dashed (blue) line, respectively, over 
our data in Figure [THI The first case is very similar to the 
case made above for the complete CNO cycle, differences 
being only due to the slightly different assumptions made 
about the compositions. 

If this scenario is correct, it would be expected that the 
observed points should lie between the two lines. Actually, 
the line representing pollution by N-poor stars does not 
reproduce the observations; on the other hand, the line 
representing N-rich stars fits the data reasonably well (al- 
though not as well the N-rich dwarfs), requiring values 
0.5 < a < 0.8 similar to those obtained in the previous 
subsection. 

The inadequacy of models with pollution by N-poor 
stars is evident when noticing that within this scheme 
we should expect to find C-poor, Na-poor stars. However, 
these stars are not observed at all (see Figure EJ- The 
inference is that in globular clusters there are no dwarfs 
polluted by RGB stars with a chemical composition typi- 
cal of field RGB stars. Within this scheme it should then 
be assumed that only stars experiencing very deep mix- 
ing polluted unevolved stars. One would then be forced 
to conclude that only a fraction of RGB stars, and only 
those in clusters, lose a great amount of mass, and that 
these very same stars do experience also very deep mix- 
ing, likely due to the same physical mechanism (rotation? 
binarity?). 

An additional problem with this scheme is that N-rich 
giants have generally no Li at their surface. We would 
then expect that Li be depleted by a factor of 2 to 5 in 
main sequence stars of globular clusters like NGC 6397, 
in contrast with observations (Bonifacio et al. 2002). 

There arc further additional concerns in a mechanism 
involving pollution by RGB stars. The lost material ends 
up polluting other stars. It cannot be a simple surface pol- 
lution: in fact, in this case there should be also noticeable 
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differences between dwarfs and subgiants (due to different 
masses of the convective envelopes) which is not observed. 
Since most of the unevolved stars observed in clusters like 
NGC 6397 and NGC 6752 are N-rich, the total amount of 
mass lost by these RGB stars should be large, about 80% 
of the cluster mass. This seems unlikely, since an RGB 
star cannot lose more than ~ 50% of its mass, the remain- 
ing being locked in the degenerate core. Another problem 
concerns the epoch when this pollution occurred. In fact, 
if the mass was lost in recent times, one would expect a 
large numbers of young stars, obviously not observed. On 
the other hand, IM-AGB stars may eject almost 80% of 
their mass (see e.g. Marigo et al. 1998), hence the mass 
requirement in this case would be met if the original initial 
mass function (IMF) of the cluster stars is not too steep, 
allowing to form many AGB stars. Evidences for a flatter 
local IMF are discussed, in this context, by e.g. D'Antona 
(2004) and Briley et al. (2001). 

6.4.2. The triple-a scenario 

As pointed out in the previous discussion, another prob- 
lem is evident from Figure^l The dilution models predict 
that stars having [0/N]-- -1.5 should have [C/N]~ -1.9, 
whereas our observations show [C/N]'^ —1.5. This sug- 
gests the presence of an additional source of ^^C, very 
likely through the triple— a process. A similar excess of 
^^C is also suggested by the ^^C/^^C isotopic ratio. Let 
us in fact assume that c is the fraction of ^^C produced 
by triple— a and (1 — c) the fraction of ^^C resulting by 
CNO-processing. Hence, we may write for ^"^C (which is 
produced only by the CNO-cycle) ^■^C = ^"^W^^ where 
Re ~ 3.5 is the equilibrium value of the isotopic ratio 
12q^13q^ We can now re-write the fraction of ^^C from 
CNO by subtracting the contribution of ^'^C and of ^^C 
by triple-a, i.e. as = (1 - c)- ^^C ^ I - c 
hence as(l — -^)x(l — c). The observed C isotopic ratio 

is then Ro = (-n^)o ~ — — ""(l-c) — ^ from which, with 

Re 

simple algebra, the fraction of ^^C due to triple— a is 

, Re 

c = 1 . 

l + Ro 

By using an observed ratio (see Table 1) of about 8, we 
derive that about 60% of the C observed is likely to come 
from triple— a burning. 

This estimate compares well with what is known from 
models of IM-AGB stars. In fact, typical estimates of the 
C abundance in NGC 6752, [C/Fe]'--^ —0.3 corresponds to 
a mass fraction of ~ 2 x 10^^ of ^^C, a value entirely con- 
sistent with model predictions for a 5 M© star of similar 
metallicity (Ventura et al. 2004). 

Up to now, our conclusions are based only on the obser- 
vation that there is too much C, with respect to the very 
large N-enhancements, and carbon isotopic ratios too high 
to be explained purely with a re-arrangement of elements 
involved in the CNO-cycle. However, other additional evi- 
dence comes from the O abundances. In the low-0, low-C 



region of Figure ^1 the existence of a sort of plateau also 
suggests that the products of triple— a are involved. This 
is not very clear in subgiants, but quite evident in dwarfs, 
whose C abundances have not been modified by the first 
dredge-up. 

7. CONCLUSIONS 

Summarizing, we can consider various mass ranges of 
likely polluters that contributed to the chemical composi- 
tion in stars presently observed: 

— metal-poor stars in the intermediate mass range, say 
1.2 < Af < 3 — 5M0, are the classical donors consid- 
ered for CH-stars (McClure and Woodsworth 1990). 
This is likely not the correct mass range for typical 
stars polluting globular clusters, since they also pro- 
duce consistent amounts of s— process elements, that 
are not seen to vary in cluster stars. 

— stars less massive than this range (the range consid- 
ered by Denissenkov and Weiss 2004) are also unlikely. 
Apart from the possible excess of ^^C discussed above, 
either (i) they release a large mass to be successively 
accreted, and this is unlikely, since the polluted stars 
would have then a large mass range, the cluster TO 
would be smeared out and we would end up with a 
cluster mostly composed of blue stragglers^; or (ii) 
they provide a small amout of polluting mass. This 
second possibility is also unpalatable because a thin 
layer accreted on the surface of a main sequence star 
would be diluted when the stars climb toward the SGB 
phase by the extension of the convective envelope. The 
evidence of large spread in C, N, O, Na abundances at 
all luminosities (see Briley, Cohen and Stetson 2002) 
strongly argues against a simple contamination of stel- 
lar surface. 

— finally the stars more massive than the above range 
are just the IM-AGB stars that, at the present, are in 
our opinion the best candidate polluters. 

Finally, we must recall that at least half or even 2/3 of 
stars observed in globular clusters seem to be heavily af- 
fected by large alterations in their chemical composition. 
In order to explain this with mass exchange from a RGB 
star onto a companion in a binary system, we would proba- 
bly end up with a huge number of low-mass X-ray binaries 
in GCs, at odds with observations. 

In summary, our findings and discussion strongly sug- 
gest that the polluters were intermediate mass AGB stars, 
and not upper RGB stars. 

On the other hand, the scenario in which IM-AGB 
stars are the primary contributors in shaping the chemi- 
cal mix of the early cluster environment still has to face 

^ Actually, we think that the mechanism proposed by 
Denissenkov and Weiss (2004) is indeed active in a minority 
of globular cluster stars belonging to binary systems, produc- 
ing the blue straggler stars with the classical McCrea (1964) 
mechanism 
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problems (Denissenkov and Herwig 2003, Denissenkov and 
Weiss 2004, Herwig 2004, Fenner et al. 2004): current mod- 
els do not seem to reproduce the required observational 
pattern. It must be however reminded that, as noticed by 
Denissenkov and Weiss, the yields computed from mod- 
els of these stars strongly depend on two poorly known 
physical inputs, namely the treatment of mass loss and 
the efficiency of convective transport (see also Ventura et 
al. 2002). Further progress in stellar modeling is strongly 
urged. 
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